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Abstract 

Neutrinoless double beta decay (0vj3j3) induced by light Majorana neutrino ex- 
change between two decaying nucleons with squark/slepton exchange inside one 
and W exchange inside the other nucleon (so-called vector-scalar exchange) gives 
£^ ■ stringent limits on R-parity violating interactions. We have extended previous work 

by including the tensor contribution to the transition rate. We discuss the improved 
limits on trilinear ^ P -MSSM couplings imposed by the current experimental limit 
on the 0vf3(3 decay half-life of 76 Ge. 
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In supersymmetric models, the new SUSY partners differ from the SM field 
content in a discrete multiplicative quantum number R-parity defined as 

Rp = (_ 1 )3B+X+25 (1) 

Here B denotes the baryon number, L the lepton number and S the spin of a 
particle leading to R p — +1 for SM particles and R p — — 1 for superpartners. 
Thus in R p conserving models superpartners can only be produced in pairs and 
the LSP is stable, leading to a natural WIMP dark matter candidate. While in 
the minimal supersymmetric extension (MSSM) of the standard model (SM) 
R-parity is assumed to be conserved, there are no theoretical reasons for R p 
conservation and several GUT [1] and Superstring [2] models require R-parity 
violation in the low energy regime. Also the reports concerning an anomaly 
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at HERA in the inelastic e + p scattering at high Q 2 and x [3] have renewed 
the interest in ^?p-SUSY (see for example [4,5]). Generally, one can add the 
following trilinear R-parity violating terms to the superpotential [1] 

Wpp = XijkLiLjEk + \ ijk LiQjD k + \ ijk U iDjD k , (2) 
where k denote generation indices, L, Q denote lepton and quark doublet 
superfields and E, U, D lepton, up- and down quark singlet superfields. Terms 
proportional to A, A violate lepton number, those proportional A violate 
baryon number. While simultaneous presence of both kinds of terms would 
lead to too fast proton decay and thus is forbidden, assuming the A" terms to 
be zero no constraints on A and A terms can be derived from proton decay. 
The search for neutrinoless double beta decay, converting a nucleus (Z, A) into 
a nucleus (Z + 2, A) under emission of two electrons, has been proven to be- 
long to the most powerful tools to search for lepton number violating physics 
beyond the SM (for a review see [6,7]). Contributions occuring through Feyn- 
man graphs involving the exchange of superpartners as well as ^?j^couplings 
A have been discussed in [8-12] and yield the most stringent bound on A m . 
Taking into account the fact that the SUSY partners of the left and right- 
handed quark states can mix with each other, also diagrams appear in which 
the neutrino-mediated double beta decay is accompanied by SUSY exchange 
in the vertices (see fig. 1). These contributions allow to constrain also combi- 
nations of couplings of higher generations A 11j A 1j1 . They have been discussed 
in [13] and more extensively in [11], where however only scalar-pseudoscalar 
currents have been taken into account, whereas the tensor contribution to the 
decay rate has been neglected. On the other hand, in a recent work [14] the 
dominance of the tensor contribution in a general framework for neutrino me- 
diated double beta decay has been proven. In the present letter we reanalyze 
SUSY-accompanied double beta decay and discuss the relative importance of 
the different nuclear matrix elements involved. 

The mixing between scalar superpartners ^ of the left and right-handed 
fermions f^R occurs due to non-diagonality of the mass matrix which can be 
written as 



m\ + mj -0A2D Z -m f (A f + ^ tan/3) 



^ -m f (A f + //tan/3) mj + mj - 0.08D Z J 

Here, f — d,s, b, e, /t, r and / are their superpartners. Dz = M§ cos2/3 with 
tan/3 = (i/aVC^i) being the ratio of vacuum expectation values of the two 
Higgs doublets, ^j LR are soft sfermion masses, Af are soft SUSY breaking 
parameters describing the strength of trilinear scalar interactions, and \x is the 
supersymmetric Higgs(ino) mass parameter. Once sfermion mixing is included, 
the current eigenstates /l,/r become superpositions of the mass eigenstates 
fi with the masses mj, and the corresponding mixing angle 9* is defined as 
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m L = 2 \ mLL + mRR T V - m RR) 2 + 4 <R) ( 4 ) 

where m 2 LR , m 2 LL , m RR denote the (1,2), (1,1), (2, 2) entries of the mass matrix 

((3))- 

Now it is straightforward to find the effective 4-fermion v — u — d — e vertex 
induced by the sfermion exchange in the diagrams presented in fig. 1. The 
corresponding effective Lagrangian, after a Fiertz rearrangement, takes the 
form 

£s f U f SY(x) = ^[\ {iMLR ~ S»L R ) ■ Ki ■ + 75)S C ) («(1 + l,)d) - 

- 2r$ )LL ■ U m • (P,(l - 75 )e, c ) (u(l + 7s )d) + (5) 
+ \v% )R r ■ Um [pi 7 "(1 + 75)e, c ) (u 7m (1 - 75 )d) + 

+ • U *m ■ ^ + ^)ei) (u <v(l + 78 )d) • 

The ^ p MSSM parameters 77 and neutrino mixing matrix C/jj are defined as 
follows 



X' X' / 1 1 

raj \ ^ ^jlfc^rafcl ^ nQ d I ' ' 



dl(fe) d 2 (k), 

A'- 1t A' i fc /sin 2 0^ cos 2 ^, 



raj _ V- ^jlfc^nlfc &m ^(fc) , uus "(ft) 1 / 7 X 



m di(fe) m <i 2 (fc) 



*, = E^gf-^,(^-^), (8) 



E^V;- ( 10 ) 



Here tj^lr denotes the contribution vanishing in the absence of /l — /r - mix- 
ing while T)tf)LL an d V(f)RR m this limit correspond to the fi and /r exchange 
contribution in fig. 1. We use the notations d(M = d,s,b and = e, u,r. 
Due to the antisymmetry of the Yukawa coupling \ n jk in nj it follows that 
VftjLR — 0- This is an essential difference between the slepton l L — l R and 
the squark — <?r contributions. The latter is not imposed to vanish at any 
combination of indexes. Since the and rj?^ RR contributions to the di- 

agram fig. 1 are helicity suppressed oc (m„) ~ C(0.5)eV, we can restrict 
ourselves to the consideration of V^lr anc ^ tfi^LR' which contributions are 
momentum enhanced oc q ~ p F ~ 100 MeV, where p F denotes the Fermi 
momentum. Here (m v ) and q denote the effective neutrino mass and momen- 
tum entering the neutrino propagator. The first line in eq. (5) corresponds to 
the scalar-pseudoscalar (S + P) contribution, the last line in (5) to the tensor 
contribution. 
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In the light neutrino case the 0u(3/3 decay rate is given by 

TyXpvW) = Goi {^ ( i)Ms + p + (-%) + %)) (M s+ P + M T )} 2 . (11) 

Here Gq\ denotes the phase space factor defined in [15], n q ) = V( q )LR anc ^ 

effective parameters are introduced as 

r K i, q) = J2^nV( l \ )L R- (12) 

n 

For the fjm summation starts from n = 2 and A n denotes the combination of 
mixing matrices corresponding to heavy or sterile neutrinos 

A* = IX-^ + t/* £U (13) 

3 

where the sum extends over heavy mass eigenstates m Vi > 10 GeV (see [11]). In 
s-wave approximation for the outgoing electrons and under some assumptions 
according to [16,11] (the s-wave approximation is expected to affect the result 
less than 10 % [14]) the matrix element is 
F^(Q) 1 

m ^ = w;( m t + - 3 m ^' < 14) 

M T = «i (^M GT > - M T >), (15) 
with (summation over nucleons a, b is suppressed) 

M ar = (0}\h R (a a a b )r:r b + \0t) " (16) 

M T > = (0}\h T/ {(a a hb)(<?bhb) ~ l(Wb)}T+T b + \0+), (17) 

,(3), 



T? ) (0)Gv (l-2m P (G w /Gy)) 
2G\Rm 

Hr and h T i are neutrino potentials defined as 



«i = 0^2 T^Z • (I 8 ) 



h R = 



7 dqq MQr ab )f 2 (q 2 ) 

7r mp J ui(u 
o v 

2 -R 2 7 J 2 Pit 



fjtf MirMto) (19) 
J u(u) + E) 



V = / dqq 2 ; [ % {q 2 3o (qr ab ) - 3^- 3l (qr ab )}. (20) 

impi tu(uj + E) r ab 

Here R denotes the nuclear radius, mp the proton mass. 
Further uj = \J<f + rn%, q = \q\, r — f/r and jk(qr) are spherical Bessel 
functions, (m + E)~ l is the energy denominator of the perturbation theory. 
The form factors F t a (0) = F^q 2 )/ f(q 2 ) and 7\ (3) (0) = T x (3) (q 2 ) / f (q 2 ) with 
f(q 2 ) = (1 + q 2 /m 2 A )~ 2 (m\ = 0.85 GeV) have been calculated in the MIT 
bag model in [17], Ga — 1.26, Gy — 1 and the strength of the induced weak 
magnetism (Gw/Gy) = ^mp" — 2~^p" * s °kt ame d by the CVC hypothesis. 
For comparison and to correct some details in [11] we first concentrate on the 
S + P part, , i.e. M.? = in eq. 11. Inserting the numerical value of the matrix 
elements M. GT > = 2.95 and Ai T ' = 0.224 [11] for the special case of 76 Ge and 
the half life limit obtained from the Heidelberg-Moscow experiment, Tf^f^ > 
1.2 • 10 25 y, [18,7] one derives r\ q < 4.5 • 10~ 8 , fj^ < 1.1 • 10~ 8 , corresponding to 
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limits of 

A„A 311 A nl3 < 9.9 • 10- 8 (^^) 3 . (21) 
311 nu _ ViooGeV^ v ; 

Here in the last equation only the term corresponding to r exchange has been 
kept. These limits differ from the ones given in [11] due to an erroneous factor 
of 2 in the definition of Ais+p m that reference as well as the improvement 
in the experimental bound. 

Keeping the tensor part in eq. 11 the half life limit of the Heidelberg-Moscow 
Experiment implies: 

O'^.lO-^)' (22) 
for supersymmetric mass parameters of order 100 GeV. For Asusy ~ 1 TeV, 
motivated by SUSY naturalness arguments, one obtains A' 112 A' 121 < 1.1 • 10~ 3 , 
A 'ii3 A 'i3i < 3-8 • 10~ 5 . These are by more than a factor of four more stringent 
than the limits obtained from the scalar pseudoscalar part considered in [11]. 
The uncertainty of the nuclear matrix elements involved can be estimated to 
be less than a factor of two. This is motivated by the fact that the 2v(3(3 
halflife of 76 Ge has been predicted correctly within a factor of 2 (the 2v(3(3 
matrix element within a factor of \/2) [19]. Since the uncertainty in both 
parts (S + P and tensor) is expected to be about the same, the improvement 
can still be considered as substantial. 

The obtained bound should be compared with the limits obtained from tree 
level K°-K° and B°-B° mixing, which yield X' il2 \' i21 < MO" 9 and A- 13 A- 31 < 
8 • 10~ 8 respectively, for m g = 100 GeV [20]. While the second generation 
is bounded by about three orders of magnitude more stringent from the K 
system, for the third generation the double beta bound is most stringent. 
Moreover, since the masses of different exchanged particles (selectrons and 
squarks) enter, the limits are complementary in some sense. 
In conclusion, we have performed a reanalysis of the SUSY-accompanied neu- 
trino exchange mode of neutrinoless double beta decay. Contrary to the pre- 
vious ansatz we included the tensor contribution to the decay rate, which has 
been shown to be the dominant contribution. This improves the limits on 
derived without the tensor contribution by a factor of four and thus 
provides the most stringent bound on A' 113 A' 131 . 
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Fig. 1. Feyman graphs for the neutrino exchange mechanism of neutrinoless double 
beta decay accompanied by (a) squark and (b) slepton exchange. 
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